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CRITICAL SURFACE OF THE 1-2 MODEL 


GEOFFREY R. GRIMMETT AND ZHONGYANG LI 

Abstract. The 1-2 model on the hexagonal lattice is a model of statistical me¬ 
chanics in which each vertex is constrained to have degree either 1 or 2. There 
are three edge-directions, and three corresponding parameters a, b, c. It is proved 
that, when a > b > c> 0, the surface given by y/a = Vb + y^is critical. The proof 
hinges upon a representation of the partition function in terms of that of a certain 
dimer model. This dimer model may be studied via the Pfaffian representation of 
Fisher, Kasteleyn, and Temper ley. It is proved, in addition, that the two-edge cor¬ 
relation function converges exponentially fast with distance when -y/a ^ Vb + y/c. 
Many of the results may be extended to periodic models. 


1. Introduction and background 

The 1-2 model on the hexagonal lattice was introduced by Schwartz and Bruck 
[32] as an intermediary in the calculation of the capacity of a constrained coding 
system. They expressed the capacity via holographic reductions (see [36]) in terms 
of the number of perfect matchings (or dimer conhgurations), and the latter may be 
studied via the Pfaffian method of Fisher, Kasteleyn, and Temperley [12, 18, 34]. 
The 1-2 model may be viewed as a model of statistical mechanics of independent 
interest, and it is related to the Ising model and the dimer model. In the current 
paper, we study the 1-2 model within this context, and we establish the exact form 
of the associated critical curve. 

A 1-2 conhguration on the hexagonal lattice El = (V, E) is a subset F of edges 
such that every vertex is incident with either one or two edges of F. There are three 
real parameters a,b,c> 0, which are associated with the three classes of edges of El. 
The weight of a configuration on a hnite region is the product over vertices v of one 
of a, b, c chosen according to the edge-configuration at v. (See Figure 2.2.) 

Through a sequence of transformations, the 1-2 model turns out to be linked to an 
enhanced Ising model, a polygon model, and a dimer model. These connections are 
pursued here, and in the linked paper [16]. The main result (Theorem 3.1) states in 

Date: 28 June 2015, revised 1 June 2016, 17 February 2017. 

2010 Mathematics Subject Classification. 82B20, 60K35, 05C70. 

Key words and phrases. 1-2 model, Ising model, dimer model, perfect matching, Kasteleyn 
matrix. 


1 



2 


GEOFFREY R. GRIMMETT AND ZHONGYANG LI 


effect that, when a > 6, c > 0, the surface given by y/a = Vb + ^/c is critical. This 
is proved by an analysis of the behaviour of the two-edge correlation function (cxecr/) 
as |e — / I —)■ oo. The model is called uniform if a = b = c= l, and thus the uniform 
model is not critical in the above sense. 

There has been major progress in recent years in the study of two-dimensional Ising 
models via rhombic tilings and discrete holomorphic observables (see, for example, [4, 
7, 8, 21]). There is a rhombic representation of the critical polygon model associated 
with the 1-2 model, and an associated discrete holomorphic function, but this is not 
explored here. 

Certain properties of the underlying hexagonal lattice are utilized heavily in this 
work, such as trivalence, planarity, and support of a action. It may be possible to 
extend many of the results of this paper to certain other graphs with such properties, 
including the Archimedean lattice (3,12^) and the square/octagon lattice (4,8^). 
Further extensions are possible to periodic models on hexagonal and other lattices. 
(See Remarks 3.2, 4.3 and Section 10.2.) 

It was shown already in [25] that a (geometric) phase transition exists for the 1-2 
model on H. An a-cluster is a connected set of vertices each having local weight a (as 
above). It was shown that there exists, a.s. with respect to any translation-invariant 
Gibbs measure, no infinite path of present edges. In contrast, for given b, c, there 
exists no inhnite a-cluster for small a, whereas such a cluster exists for large a. The 
a.s. uniqueness of inhnite ‘homogeneous’ clusters was proved in [27]. 

This paper is concentrated on the 1-2 model and its dimer representation. A 
related representation involves the polygon model on H, and the phase transition of 
the latter model is the subject of the linked paper [16]. The polygon representation 
is related to the high temperature expansion of the Ising model, and results in an 
inhomogeneous model that may regarded as an extension of the 0{n) model with 
n = 1; see [11] for a recent reference to the 0{n) model. 

The structure of the current work is as follows. The precise formulation of the 1-2 
model appears in Section 2, and the main theorem (Theorem 3.1) is presented in 
Section 3. 

The 1-2 model is coupled with an Ising model in Section 4, in a manner not 
dissimilar to the Edwards-Sokal coupling of the random-cluster model (see [15, Sect. 
1.4]). It may be transformed into a dimer model (see [25]) as described in Section 
5. In Section 6, we gather some conclusions about inhnite-volume free energy and 
inhnite-volume measures that are new for the 1-2 model. Theorem 3.1 is proved in 
Sections 7-8 by an analysis using Pfafhans, and further in Sections 4.5 and 9. Section 
10 is devoted to extensions of the above results to periodic 1-2 and Ising models to 
which the Kac-Ward approach of [29] does not appear to apply. 
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2. The 1-2 model 

Let G = iV,E) be a finite graph. A 1-2 configuration on G is a subset FOE 
such that every n G iL is incident to either one or two members of F. The subset 
F may be expressed as a vector in the space S = {—where —1 represents 
an absent edge and -|-1 a present edge. Thus the space of 1-2 configurations may be 
viewed as the subset of S containing all vectors a such that 

^<€{ 1 , 2 }, veV, 

eBv 

where 

(2.1) a'(e) = i(l + a(e)). 

(In Section 4.2, we will write E® for E, in order to distinguish it from a space of 
vertex-spins to be denoted E'".) 



Figure 2.1. An embedding of the hexagonal lattice. Horizontal edges 
are said to be of type a, NW edges of type b, and NE edges of type c. 

Suppose now that G is a finite part of the hexagonal lattice H, suitably embedded 
in see Figure 2.1. The embedding is such that each edge may be viewed as one of: 
horizontal, NW, or NE. (Later we shall consider a finite box with toroidal boundary 
conditions.) Let a, 6, c > 0 be such that (a, b, c) ^ (0,0, 0), and associate these three 
parameters with the edges as indicated in the figure. For a G E and v G V, let a\^ 
be the sub-configuration of a on the three edges incident to v. There are 2^ = 8 pos¬ 
sible local configurations, which we encode as words of length three in the alphabet 
with letters {0,1}. That is, for v E V, we observe the states 

where Cy^a, ^v,b) ^v,c are the edges of type a, b, c (respectively) incident to v. The 
corresponding signature Sy is the word cT'(e^,c)<T'(e^,6)cT'(e^,a) of length 3, where a' is 
given in (2.1). That is, the signature of v is given as in Figure 2.2, together with the 
local weight tc(cr|^) associated with each of the eight possible signatures. 

The hexagonal lattice El is, of course, bipartite, and we colour the two vertex- 
classes black and white. The upper diagrams of Figure 2.2 are for black vertices, and 
the lower for white vertices. 
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Figure 2.2. The eight possible local conhgurations a\v at a vertex v 
in the two cases of black and white vertices. The signature of each is 
given, and also the local weight tc(cr|^) associated with each instance. 


To the vector a G S, we assign the weight 

(2.2) w(cr|^). 

vGV 

These weights give rise to the partition function 

(2.3) Z = 

(tGS 


which leads in turn to the probability measure 

(2.4) /i(a) = ^w{a), a G S. 

It is easily seen that the measure fi is invariant under the mapping (a, b, c) h-)■ 
{ka, kb, kc) with k > 0. It is therefore natural to re-parametrize the 1-2 model 
by 


(2.5) 


(a',6',c') 


(,a,b,c) 

\\{a,b,c)\\ 2 ' 


We will work mostly with a hnite subgraph of El subject to toroidal boundary 
conditions. Let n > 1, and let ri, T 2 be the two shifts of H, illustrated in Figure 2.3, 
that map an elementary hexagon to the next hexagon in the given directions. The 
pair (ti, T 2 ) generates a action on H, and we write EI^ for the quotient graph of El 
under the subgroup of generated by r” and r^. The resulting EI„ is illustrated in 
Figure 2.3, and may be viewed as a hnite subgraph of El subject to toroidal boundary 
conditions. 
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Figure 2.3. The graph EI„ is an n x n ‘diamond’ wrapped onto a 
torus, as illustrated here with n = 4. 

Our purpose in this paper is to study the 1-2 measure (2.4) on ]HI„ in the inhnite- 
volume limit as n —)■ cxd, and to identify its critical surface. As an indicator of phase 
transition, we shall use the two-point function (aea/)„, where e, / are two edges and 
{■)n denotes expectation. 

We do not explore in detail the nature and multiplicity of inhnite-volume measures 
in this paper. There are certain complexities in such issues arising from the absence 
of a correlation inequality, and some partial results along these lines may be found 
in [25, Thm 0.1]. These results are developed in Section 6, where the main result of 
current value is the existence of the inhnite-volume limit of the toroidal 1-2 measure, 
see Theorem 6.2. 


3. Main results 

Consider the 1-2 model on EI^ with parameters a,b,c > 0. We write e = {x, y) 
for the edge e with endpoints x, y, and we use {X)n to denote expectation of the 
random variable X with respect to the probability measure of (2.4) on ]HI„. We 
shall make use of a measure of distance \e — f\ between e and /, and it is largely 
immaterial which measure we take. For dehniteness, consider El embedded in in 
the manner of Figure 2.3, with unit edge-lengths, and let |e — f\ be the Euclidean 
distance between their midpoints. 

We shall sometimes require the following geometric condition on two NW edges 
e,/ e E; 


(3.1) 


there exists a path tt = 7r(e, /) of EI„ from e to / 
using only horizontal and NW half-edges. 
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Theorem 3.1. Let a,b,c> 0, and e, / G E. 

(a) The limit {(Te^f) = limn-)-oo(o'eO'/)n exists. 

(b) Subcritical case. Let a > 6 > 0 and ^/a — y/b < yfc < yfa + a/6. There exists 
a{a, b,c) > 0 such that 

(3.2) |{<T.<T,)| < e,/eE. 

(c) Supercritical case. Let a>b>0, and let e, f be NW edges satisfying (3.1). 
For almost every c > 0 satisfying either ^/a > y/b + ^/c or ^/c > ^/a + y/b, 
we have that lim|e-/|^.oo{o'eO'/)^ exists and is non-zero. The convergence is 
exponentially fast in the distance \e — f\. 

The two-edge function {(TeCTf) behaves (when a > 6 > 0) in a qualitatively different 
manner depending on whether or not ^/a — y/b < a/c < ^/a+y/b. Here is a motivation 
for condition (3.1). Consider the ‘ground states’ when either c = 0ora = 6 = 0. By 
examination of the different cases in Figure 2.2, we may see, subject to (3.1), that 

(3.3) {o'eO'f) = 1 if either a, 6 > 0, c = 0, or a = 6 = 0, c > 0. 

The result of part (c) will follow from this by an argument using analyticity (and, 
moreover, the set of c at which the conclusion of (c) fails is a union of isolated points). 
Part (c) holds with e, / assumed to be horizontal rather than NW. 

Theorem 3.1 is not of itself a complete picture of the location of critical phenomena 
of the 1-2 model, since the conditions on the parameters in part (c) are allied to 
the direction of the vector from e to /. (The direction NW is privileged in the 
above theorem. Similar results hold for the other two lattice directions with suitable 
permutations of the parameters.) We have not ruled out the theoretical possibility 
of further critical surfaces in the parameter-space [0, cx))^. 

Remark 3.2. The guickest proof of Theorem 3.1(h), the subcritical case, (given in 
Section 4-5) is based on a result of [29] that imposes a condition on the parameters 
of edges incident to a vertex v, uniformly in v. This condition is satisfied in the 
current setting (see Section 4-5). In the more general setting of certain periodic but 
non-constant families of parameters, or possibly of the 1-2 model on other graphs 
such as the sguare/octagon lattice, much of Theorem 3.1 remains true, but the con¬ 
dition of [29] does not generally hold. In order to overcome this lacuna for more 
general systems, we present a further proof of Theorem 3.1(h) in Section 9 (in the 
more general form of Theorem 9.1) using the dimer-related technigues of the proof of 
Theorem 3.1. Such results may be extended in part to more general periodic settings, 
see Section 10.2. 

The proof of Theorem 3.1 utilizes a sequence of transformations between the 1-2 
model and the Ising and dimer models, as described in the forthcoming sections. 
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Theorem 3.1(b) is proved in Section 4.5. Most of the remaining proof is found in 
Section 7, with the exponential rate of part (c) proved in Section 8. The last is 
proved via a general result concerning the convergence rate of the determinants of 
large truncated block Toeplitz matrices to their limit when the symbol is a smooth 
matrix-valued function on the unit circle. 


4. Spin representations of the 1-2 model 

Two spin representations of the 1-2 model are presented here. In the first, the 1-2 
partition function is rewritten in terms of edge-spins. The second is reminiscent of 
the random-cluster representation of the Potts model. A further set of spin-variables 
are introduced at the vertices of the graph, together with an Ising-type partition 
function. 


4.1. The 1-2 model as a spin system. Let EI^ = {Vn, En) be the quotient hexago¬ 
nal lattice embedded in the torus in the manner of Figure 2.3. Let = {—1, -t-1}®", 
where —1 (respectively, -|-1) represents an absent edge (respectively, present edge). 

For cr G and v G 14, let cr^,6, <t„,c denote the spins on the incident a- 
edge, 5-edge, c-edge of v. Two partition functions Z, Z' generate the same measure 
whenever they differ only in a multiplicative factor (that is, their weight functions 
satisfy w(a) = cw'(a) for some c 4 0 and all a G S), in which case we write Z = Z'. 
We represent the 1-2 model as a spin system as follows. 

Proposition 4.1. Let a,b,c > 0 such that {a,b,c) 4 (0,0,0). The 1-2 model with 
parameters a, b, c on ]HI„ has partition function Zn satisfying Z^ = Z'^ where 


(4.1) 


n (f T T E (J y^ “t“ , 


CtGSti vGVn 


and 

(4.2) 


A 


a — b — c 
a + b + c 


B 


b — a — c c — a — b 

1 7 * 

a-l-o-fc a-l-o-l-c 


Proof. By examination of (4.1), we see that a vertex with local configuration labelled 
a in Figure 2.2 has weight 


1+A-B-C 


4a 

a -|- 5 -|- c’ 


with similar expressions for vertices with the other possible signatures. This is in 
agreement with (2.2)-(2.3), and the claim follows. □ 
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4.2. Coupled Ising representation. Let AM.n = (A14, AEn) be the graph derived 
from EI„ = (VnyEn) by adding a vertex at the midpoint of each edge in En. Let 
MEn = {Me : e G E^} be the set of such midpoints, and AVn = 14 U MEn- The 
edges AEn are precisely the half-edges of En, each being of the form {v, Me) for some 
n G 14 and incident edge e E En- 

We introduce an Ising-type model on the graph 4EI„. The marginal of the model 
on midpoints MEn is a 1-2 model, and the marginal on I 4 is an Ising model. This 
enhanced Ising model is reminiscent of the coupling of the Potts and random-cluster 
measures, see [15, Sect. 1.4]. It is constructed initially via a weight function on 
conhguration space, and via the associated partition function. The weights may 
be complex-valued, and thus there does not always exist an associated probability 
measure. 

The better to distinguish between I 4 and MEn, we set S® = {—1,-1-1}^'®’* as 
before, and = { — 1, -fl}^". An edge e E En is identihed with the element of MEn 
at its centre. A spin-vector is a pair a^) G E® x S'" with cr® = (cr^,s : v E 14, s = 
a,b,e) and ex'' = : v E 14), to which we allocate the (possibly negative, or even 

complex) weight 

(4.3) n (f T enCTydyn) i^i- T ^b^v^v,b) A ^c^v^v^c): 

V&Vn 

where e^, e^, ec E C are constants associated with horizontal, NW, and NE edges, 
respectively, and dy^a, o'yfi, dy^c denote the spins on midpoints of the corresponding 
edges incident to n G 14- If u and v are endpoints of the same edge {u, v) of EI„, then 
o'u,a = o'y^a- In (4.3), each factor 1 -|- egdydy^s (s = a, b, c) corresponds to a half-edge 
of EI„. Recalling that 

(4.4) _|_ (j^(j 2 tanhx) coslix, a: G M, did 2 = ±1, 

the above spin system is a ferromagnetic Ising model on AEI„ when Cq, €(,, €c E (0,1). 

4.3. Marginal on the midpoints MEn- The partition function of (4.3) is 

(4.5) Zn{I) : = E E n (1 -F eadydy^a){i- + ebdydy^b){l + e^dydy^^. 

0-®eSe O-'^eS^ vEVn 

(The notation Zn{I) is chosen for consistency with the polygon partition function 
Zn{P) used in this article and imported from [16].) The product, when expanded, 
is a sum of monomials in which each dy has a power between 0 and 3. On summing 
over a'", only terms with even powers of the site-spins dy survive, and furthermore 
= 1, so that 

Zn{I) = 2l^"l E n (1 “1“ ^b^c^v,b^v,c ^a^c^v,a^v,c ^a^b^v^a^v,b^ • 

VGVn 
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Let a, 6, c > 0 be such that ABC ^ 0 where A, B, C are given by (4.2), and let 



By (4.1), 

(4.7) Z^{I) = 

whence the marginal model of (4.3) on the midpoints of edges of EI^, subject to (4.6), 
is simply the 1-2 model with parameters a, b, c. 


4.4. Marginal on the vertices Vn- This time we perform the sum over a® in (4.5). 
Let g = {u, v) G En be an edge with weight eg. We have 

(4.8) (1 + egCr^cTg) = 2 [l + ejcr^a^) , 

(Tg=±l 

(4.9) ^ ag{l + egau<Jg)(l + egayag) =2eg{au + <Jy). 

crg=±l 

By (4.5) and (4.8), 

(4.10) Z„(/) = 2l®"l H {1 + elauay). 

o-''eS^ g=(u,v)£E„ 


By (4.4), this is the partition function of an Ising model on ]HI„ with (possibly com¬ 
plex) weights. 

Let e = {u,v), f = {x,y) be distinct edges in En- Motivated by Section 4.3 and 
the discussion of the two-edge correlation (creCr/)„ of the 1-2 model, we dehne 


(4.11) a(e,/) 


1 


E E (1 T eaCT^jCT^^a) (1 T CfeCTtiCr^i^b) (1-f-ecCTt,CT^,^c) • 

O-eeSe veVn 


By (4.9)-(4.10), this equals 


Znil) 


2\E„\ ^ ^e{<^u + 0'v)^fi0'x + (^y) 

JtC^+ (1+ '‘A-A 5.(hhs 




^ u;(ct''). 




<7''esv 


(4.12) 
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where 


w{a^) 




g={u,v)^En 

^eis^u fij^x H“ ^y) 

(1 + elau(Jv) (l + eja^ay) ’ 


e = {u,v), f = {x,y). 


We interpret Dej^cr'') as 0 when its denominator is 0. Since cxi + <72 = 0 when 
o'iO '2 = —1, we may write 


(4.13) 




“1“ d" ^y') 

{l + el){l + e}) ' 


If the weights w{a'') are positive (which they are not in general), the ratio on the 
right side of (4.12) may be interpreted as an expectation. This observation will be 
nsed in Section 4.5. 

By inspection of (4.10), if = ±1 for some g G {a, b, c}, then zero mass is placed 
on conhgnrations a for which there exists an edge {u,v) of type g with cr„cr^ = =f1 - 
We tnrn to the special case of (4.6) and (4.2) with ABC ^ 0. Then 


(4.14) 


-1 if and only if = 6^ + c^ 
1 if and only if 6c = 0, 


and similarly for eb, Cc. Note in this case that a{e, f) = {ae<Jf)n, the two-edge 
function for the associated 1-2 model. 


4.5. Proof of Theorem 3.1(b). Let e = {u,v), f = {x,y) be distinct edges in En 
snch that u, x are white and v, y are black. By (4.11)-(4.13), 

(4-15) = {Dej{a^))\, 

where (■)Jj denotes expectation in the Ising model of (4.10). Recall that this Ising 
model may have complex weights. 

By [29, Cor. 2.5] and known results for the Kac-Ward operator (see [6, 17, 28]), 
we have that {o'e^'f) ■= liiiin-)-oo(<^eO'/)n —t 0 exponentially fast as \e — f\ —)■ oo, so 
long as the three acute angles with tangents |e^|, g = a,b,c, have sum 9 satisfying 
9 < \n. 

It suffices to assume that a > 6, c > 0 and i/a < a/6 -|- -/c. Suppose hrst that, in 
addition. 


(4.16) 


a / 6 -|- c. 
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Let A, B, C be given by (4.2) and (4.6), so that A G (—1,1) \ {0} and B^C < 0. 
Note that the eg of (4.6) are purely imaginary if ^4 < 0, and real otherwise. Now, 

(4.17) tang=|Ai^C| ^_ 

which is hnite under (4.16) and strictly positive if 

(4.18) A^ + B‘^ + C^ < 1. 

Using (4.2), it is a short calculation to see that (4.18) holds if — 2ah — 

2hc — 2ca < 0, which is indeed valid when < y/b -V yjc. (See also the proof of 
Proposition 5.1.) This establishes (3.2) subject to (4.16). 

Suppose hnally that a = 6 + c, so that A = 0 and B,C < 0, B + C = —1. It 
is useful to represent the 1-2 model as a polygon model, via its high-temperature 
expansion. As explained in [16], the two-edge function satishes 


{(Je<Jf)n 


Zn{P) ’ 


where Zn,e-^f and Zn{P) are given at [16, eqns (2.3), (2.7)] with 


2l, P} 


For a polygon conhguration vr (that is, a set of edges such that every vertex has 
even degree), a vertex of ]HI„ is said to be of type ab if it is incident to two edges with 
types a and b (and similarly for ac and be). Since each vertex in the polygon model 
has even degree, and A = 0, no vertex of EI^ has type be. Therefore, any polygon 
conhguration with non-zero weight in Zn{P) is a disjoint union of cycles comprising 
ac-type and ab-tjpe vertices. The vertices on such a cycle form consecutive pairs with 
the same type, and each such pair contributes weight either B^ or It follows that 
Zn{P) is a sum of positive weights. 

Suppose hrst that e and / are a-type (horizontal) edges. Let tt' be a path between 
the midpoints of e and / that contributes a non-zero weight to Zn,e<^f, and let h be 
the number of its a-type edges (with each a-type half-edge contributing 7). Then vr' 
contains exactly 2h vertices of ]HI„, which appear in consecutive pairs with the same 
type (either ab or ac). The product of the weights of the vertices of vr' is B‘^'" 

(> 0), where v is the number of consecutive pairs with type ac. We denote by T(h, v) 
the set of all such vr'. 
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Since the removal of vr' gives a configuration contributing to Zn{P), and in addition 


D:= B^ + C^ < 1, 




where H = inf{h : T{h,v) 7^ 0 for some v}. Since H > c\e — f\ for some c > 0, the 
claim follows for horizontal e, /. 

If either e or / is not horizontal, an extra term appears at one or both of the ends 
of tt', and such a term contributes a factor bounded by | max{i?, C}| < 1. 

Remark 4.2. The conclusion (3.2) of Theorem 3.1(h) may he proved as follows 
subject to the more restrictive condition of . Under this condition, we have 

that A,B,C < 0. The graph EI^ is bipartite with vertex-classes coloured black and 
white (see the discussion around Figure 2.2). We now reverse the signs of the spins 
of black vertices, thereby obtaining a ferromagnetic Ising model. It is easily checked 
that this is a high-temperature model (as in (4.17)^), and it follows that (o'eO'f) —)■ 0 
exponentially fast as \e — f\ ^ 00 . 

Remark 4.3. A further proof of parts of Theorem 3.1 is presented in Section 9. 
This proof is based on ‘dimer’ rather than ‘Ising’ methods, and may he extended to 
periodic 1-2 models which appear to be currently beyond the Kac-Ward technigues of 


[29]. See Section 10. 


5. Dimer representation of the 1-2 model 


5.1. The decorated dimer model. Let = (fn,A 5 be the decorated 

toroidal graph derived from EI„ and illustrated on the right of Figure 5.1. It was 
shown in [25] that there is a correspondence between 1-2 conhgurations on EI,^ and 
dimer conhgurations on EI„ a- This correspondence is summarized in the hgure 
caption, and a more detailed description follows. 

Let cr be a 1-2 conhguration on EI„, and let n G 14 (4 14 ,a)- The vertex v has 
three incident edges in EI„^a, which are bisectors of the three angles of EI„ at v. Such 
a bisector edge is present in the dimer conhguration on EI„ a if and only if the two 
edges of the corresponding angle have the same cr-state, that is, either both or neither 
are present. The states of the bisector edges determine the dimer conhguration on 
the entire EI„_a- Note that the 1-2 conhgurations a and —a generate the same dimer 
conhguration, denoted D„. 
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Figure 5.1. Part of a 1-2 configuration on EI„, and the corresponding 
dimer (sub)configuration on EI„ a- When two edges with a common 
vertex of EI„ have the same state in the 1-2 model, the corresponding 
‘bisector edge’ is present in the dimer conhguration. The states of 
the bisector edges determine the dimer conhguration on the rest of 
EI„_a- The edges of EI„_a are allocated weights consistently with the 
1-2 weights of Figure 2.2. The central lozenge of the right-hand hgure 
is expanded in Figure 5.3. 


To the edges of Hji a we allocate weights as follows: edge e = (i,j) is allocated 
weight Wij where 


if e is a horizontal bisector edge, 
if e is a NW bisector edge, 
if e is a NE bisector edge, 
otherwise. 

The weight of a dimer conhguration is the product of the weights of present edges. 

To each 1-2 conhguration a on EI„, there corresponds thus a unique dimer con¬ 
hguration on EI„_a- The converse is more complicated, and we preface the following 
discussion with the introduction of the planar graph derived from EI„ by a process 
of ‘unwrapping’ the torus. 

Let EI(^ be the planar graph obtained from EI„ by cutting through the two homology 
cycles ■jx and 'jy of the torus, as illustrated in Figure 5.2. That is, may be viewed 
as the set of edges that intersect the region marked in Figure 5.2 (in which n = 4 
and the central edge is labelled {u,v)). We may consider EI(^ as a ‘partial-graph’ 
EI'„ = (Vn, En, Hn), wliere W is the vertex set, is the ‘internal’ edge set, and Hn 
is the set of half-edges having one endpoint in W and one outside W- We write H^. 
and (respectively, Hy, Hy) for the sets of half-edges that cross the upper left and 


(5.1) 






= 
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lower right sides (respectively, upper right and lower left sides) of the diamond of 
Figure 5.2. Let Hu = U for u = x,y. 



Figure 5.2. The n x n ‘diamond’ of H, with n = 4. The region 
comprises all edges and half-edges that intersect the larger diamond. 
The annulus between the given boundaries comprises a cycle (drawn 
above in bold), and two further edges incident with the Wi. 


A 1-2 conhguration on is a subset of edges and half-edges such that, for u G 14, 
the total number of edges and half-edges that are incident to v is either 1 or 2. 
It is explained in [25, p. 4] that dimer conhgurations on EI„ a are in one-to-two 
correspondence to 1-2 conhgurations on EI(^ satisfying any of the following (pairwise 
exclusive) conditions: 

(ss) for e G H^UEy, the two corresponding half-edges G H^DHy, G 
have the same state (either both are present or neither is present); 

(os) for e E H^, the two corresponding half-edges G H^, G have the oppo¬ 
site states (exactly one of them is present); for e G Hy, the two corresponding 
half-edges E Hy, E Hy have the same state; 

(so) for e E Hx, the two corresponding half-edges G H]., G H^ have the same 
state; for e G Hy, the two corresponding half-edges G Hy, E Hy have 
the opposite states; 

(oo) for e G the two corresponding half-edges G HlVJHy, E H^UEy 

have the opposite states. 

We refer to the above as the mixed boundary condition on EI 4 

The above mixed boundary condition is more permissive than the periodic con¬ 
dition that gives rise to 1-2 conhgurations on the toroidal graph EI„, although the 
difference turns out to be invisible in the inhnite-volume limit (see Theorem 6.2). 
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5.2. The spectral curve of the dimer model. We turn now to the spectral curve 
of the above weighted dimer model on EI„ a- The reader is referred to [26] for relevant 
background, and to [25, Sect. 3] for further details of the following summary. 

The fundamental domain of EI^ a is the central lozenge of Figure 5.1, as expanded 
in Figure 5.3. The edges of EI„ a are oriented as in the latter hgure. It is easily checked 
that this orientation is ‘clockwise odd’, in the sense that any face of ]HI„^a, when 
traversed clockwise, contains an odd number of edges oriented in the corresponding 
direction. The fundamental domain has 16 vertices, and its weighted adjacency 
matrix (or ‘Kasteleyn matrix’) is the 16 x 16 matrix B = {bij) with 

{ Wij if {i,j) is oriented from i to j, 

—Wij if {i,j) is oriented from j to i, 

0 if there is no edge between i and j, 

where Wij is given by (5.1). From B we obtain a modified adjacency (or ‘Kasteleyn’) 
matrix B(x,tc) as follows. 



Figure 5.3. A single fundamental domain of the decorated graph 
Eln.A obtained from the central lozenge of Figure 5.1. See that hgure 
for an illustration of the relationship between this fundamental domain 
and the original hexagonal lattice H. Note the homology cycles 
of the torus, and also the two weight-c edges crossed by the central 
dashed line. 

We may consider the graph of Figure 5.3 as being embedded in a torus, that is, we 
identify the upper left boundary and the lower right boundary, and also the upper 
right boundary and the lower left boundary, as illustrated in the hgure by dashed 
lines. 
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Let w, z G C be non-zero. We orient each of the four boundaries of Figure 5.3 
(denoted by dashed lines) from their lower endpoint to their upper endpoint. The 
‘left’ and ‘right’ of an oriented portion of a boundary are as viewed by a person 
traversing in the given direction. 

Each edge (m, v) crossing a boundary corresponds to two entries in the weighted 
adjacency matrix, indexed {u,v) and {v,u). If the edge starting from u and ending 
at V crosses an upper-left/lower-right boundary from left to right (respectively, from 
right to left), we modify the adjacency matrix by multiplying the entry {u,v) by z 
(respectively, z~^). If the edge starting from u and ending at v crosses an upper- 
right/lower-left boundary from left to right (respectively, from right to left), in the 
modihed adjacency matrix, we multiply the entry by w (respectively, w~^). We 
modify the entry [v, u) in the same way. The ensuing matrix is denoted B(z, w), for 
a dehnitive expression of which, the reader is referred to [25, Sect. 3]. 

The characteristic polynomial is given (using Mathematica or otherwise) by 

(5.2) P(z, w) = det B( 2 :, w) = /(a, b, c; w, z), 


where 

/(a, 5, c;w,z) 


+ c^ + -|- 6a^c^ -|- 6b‘^c^ — 2ab ^z H—^ (a^ + b'^ — c^) 

-2ac(w + -) (a^ + c^ - b^) - 2bc (- +-) (b^ + c^ - a^) . 

\ w J ^ ' \w z / ^ 


The spectral curve is the zero locus of the characteristic polynomial, that is, the 
set of roots of P{z,w) = 0. It is proved in [25, Lemma 3.2] that the intersection 
of P{z,w) = 0 with the unit torus is either empty or a single real point (1,1). 
Moreover, in the situation when P(l, 1) = 0, the zero (1,1) has multiplicity 2. It 
will be important later to identify the conditions under which P(l, 1) = 0. 


Proposition 5.1. Let a, 6 > 0 and c > 0. 

(a) If any of the following hold, 

(i) ■s/a = \/bP a/c, (ii) Vb= y/c+ -y/a, (hi) yfaP Vb, 

the curve P{z, w) = 0 intersects the unit torus = {{z, w) : \z\ = 1, |tc| = 1} 
at the unique point (1,1). 

(b) If none of (i)-(iii) hold, the curve does not intersect the unit torus. 

Proof. The intersection of P{z, w) = 0 with can only be either empty or a single 
point (1,1), by [25, Lemma 3.2]. Moreover, since 

(5.3) /(a, b, c; 1,1) = (a^ -|- 5^ -|- — 2ab — 2bc — 2ca)^, 

we have that /(a, b, c; 1,1) = 0 if and only if y/a ± a/O ± ^/c = 0. 


□ 
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We note for future use that 


(5.4) P(l, 1) — /(a, 6, c; 1,1) — ^ + C'^ — l){a + b + c)^] , 


where A, B, C are as in (4.2). 


6. Infinite-VOLUME limits 


This paper is directed primarily at the asymptotic behaviour of the two-edge cor¬ 
relation function of the 1-2 model, rather than at the existence and multiplicity of 
inhnite-volume measures. Partial results in the latter direction are reported in this 
section. In Section 6.1, the weak limit of the toroidal 1-2 measure is proved via a 
relationship with the dimer model on a decorated graph. In Section 6.2 we prove 
the non-uniqueness of Gibbs measures for the ‘low temperature’ 1-2 model. The 
existence of the inhnite-volume free energy is proved in Section 6.3. 

6.1. Toroidal limit measure. The 1-2 model may be studied via the dimer repre¬ 
sentation of Section 5. The dimer convergence theorem of [25] is as follows. 

Theorem 6.1. [25, Prop. 3.3] Consider the dimer measure 6n,A on with pa¬ 
rameters a,b,c> 0. The limit measure 6 a ■= hm„_,.oo 5n,A exists and is translation- 
invariant and ergodic. 

Let (respectively, fin) be the 1-2 probability measure on (respectively, on 
the toroidal EI^) with parameters a, b, c and mixed boundary condition. By the 
results of [25] and the invariance of /i™“ under sign changes. 



( 6 . 1 ) 


where is the dimer conhguration on EI^^a corresponding to the 1-2 conhgura- 
tion a on EI(^. Since the topology of weak convergence may be given in terms of 
hnite-dimensional cylinder events, the weak convergence 6n,A 5a entails the weak 
convergence of /i™“ to some probability measure on El. By Theorem 6.1, 
is translation-invariant. It is noted at [25, p. 17] that the ergodicity of 5 a does not 
imply that of and indeed there exist parameter values for which /i™“ is not 
ergodic, by the result of [25, Thm 4.9]. 

Theorem 6.2. Let a,b,c > 0. The limit fi^o ■= hm„_,.oo/in exists and satisfies 
Too = /i““- In particular, for edges e, f ofM, the limit 


( 6 . 2 ) 


(fTfifT/) — lim (fT.fT/)„ 


exists. 
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Proof. Let Qn,A be the sample space of the dimer model on EIji^a- Let be the 
probability measure of the dimer model on ]HI„ a on the subspace of conhgu- 
rations with the property that, along each of the two zigzag paths of El that are 
neighbouring and parallel to and there are an even number of present bisector 
edges. 

As explained above (see also [25]), elements of fin ,a correspond to 1-2 model 
conhgurations on with the mixed boundary condition, and of fl^/^ to 1-2 model 
conhgurations on the toroidal graph ]HI„. We show next that 

(6.3) <5:)a ^ <5a, 

where (5 a := hmn^.oo ^n,A is given in Theorem 6.1. 

Let Z^^A (respectively, be the partition function of fin ,a (respectively, 

and let Kn{z, w) be the modihed Kasteleyn matrix of ]HIn,A (see [25] and Section 5.2). 
As explained in [32, Sect. 4B], for z,w E { — 1,1}, Pf Kn{z, w) is a linear combination 
of partition functions of dimer conhgurations of four different classes, depending on 
the parity of the present edges along the two zigag paths winding around the torus. 
In particular, by [32, Table 1, Sect. 4B], when n is even, 

Pf iL„(l, 1) + Pf iLn(-l, 1) + Pf lPn(l, -1) + Pf lPn(-l, "1) • 

Let Cj = (wj, fj), 1 < i < A;, be edges of Eln^A, and let M(ei,..., e^) be the event that 
every Cj is occupied by a dimer. Let tCj > 0 be the edge weight of e*. Then 


^ev 1 

^n,A 4 


(5r(M(ei,...,efc)) 

k 

i=l 

where Kn is the submatrix of Kn obtained by removing rows and columns indexed 
by Ml, ui,..., Uk, Vk- As in [3, Thm 4], 

(6.4) 


Pf iP^(l, 1) + Pf iLn(-l, 1) + Pf iPn(l, -1) + Pf iPn(-l, -1) 
Pf iP„(l, 1) + Pf A:.(-1, 1) + Pf iPn(l, -1) + Pf lPn(-l, -1) 


(5n,A(Af(ei,... ,efe)) 

-Pf iL^(l, 1) + Pf iPn(-i, 1) + Pf iPn(i, -1) + Pf iPn(-i, -1) 

-Pf A:„(1, 1) + Pf iL.(-l, 1) + Pf iPn(l, -1) + Pf iPn(-l, -1) ■ 

As in the proof of [3, Thm 6], ..., e^)) and 5n,A(Af(ei,..., e^)) converge 

as n —>■ cx) to the same complex integral. Since the events M{ei,, Ck) generate the 
product (T-held, we deduce (6.3). 

Finally, we deduce the claim of the theorem. An even (respectively, odd) corre¬ 
lation function is an expectation of the form {a a), with a a = YleeA^e where A is 
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a finite set of edges of H with even (respectively, odd) cardinality. In order that 
fin —/i™“, it suffices that the correlation functions of and have the same 
limit. By invariance under sign change, the odd correlation functions equal 0. 

The relationship between a 1-2 measure /i and the corresponding dimer measure 
6 is as follows. Let Ci,..., Cfc be bisector edges of EI^ a, and let S{ei ,..., e^) be the 
event that every e, separates two edges of EI^ with the same 1-2 state. Using the 
correspondence between 1-2 and dimer conhgurations, 

fi{S{ei ,..., Cfe)) = h(M(ei,..., Cfc)). 

Let k > 1, let Ci, 62 ,..., e 2 k be distinct edges of EI„, and write cij = erg.. Then 

(6.5) (di ■ ■ ■ a2k)^l = 1 - 2p(cri ■ ■ ■ a2k = -!)• 

For i = 2,3,..., 2k, let vr* be a self-avoiding path between the midpoints of ei and e* 
comprising edges of EI„ and two half-edges, and let Ai be the event that the number 
of absent bisector edges encountered along vij is odd. As we move along Tii in the 
1-2 model, the edge-state changes at a given vertex if and only if the corresponding 
bisector edge is absent. Therefore, cricxi = —1 if and only if Ai occurs, so that 

(6.6) ((TiCTi)^ = fi{Ai) - fi{Ai). 

Let A be the event that the set I = {i : Ai occurs} has odd cardinality. Since 
I = {i ■. ai ^ ai}, we have that 

(6.7) ij,(ai---a2k =-1) = 5{A). 

We return to the measures fin and By (6.3)-(6.7), the even correlation 

functions of fin and are convergent as n —)■ 00 , with equal limits. It follows that 
fin fioo where fi^ = □ 

6.2. Non-uniqueness of Gibbs measures. We show the existence of at least two 
Gibbs measures (that is, ‘phase coexistence’) for the ‘low temperature’ 1-2 model 
on EI. Let S be the set of 1-2 conhgurations on the inhnite lattice EI, and let 
Q = G{a,b,c) be the set of probability measures on S that satisfy the appropriate 
DLR condition. (We omit the details of DLR measures here, instead referring the 
reader to the related discussions of [ 2 , Sect. 2.3] and [15, Sect. 4.4].) Since S is 
compact, by Prohorov’s theorem [ 1 , Sect. 1.5], every sequence of probability measures 
on S has a convergent subsequence. It may be shown that any weak limit of hnite- 
volume 1-2 measures lies in Q, and hence Q ^ 0. 

Theorem 6.3. Let a > b > 0. For almost every c satisfying either Q < ^/c < ^/a—\/b 
or yA > ^/a -|- a/5, we have that \Q\ > 2. 

Proof. This proof is inspired by that of [2, Thm 6.2], and it makes use of Theorem 
3.1, the proof of which has been deferred to Section 7. Let e be a given horizontal 
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edge of EI„, and let fm be an edge satisfying (3.1) snch that 7 r{e, fm) has length m. 
By Theorem 3.1(b) and translation invariance, for almost every c satisfying the given 
inequalities, there exists a > 0 such that 

(6.8) lim {aeCTfJ^ = 

m^oo 

where {ceCf) is the limiting two-edge correlation as n —>■ oo (see Theorem 6.2). It 
suffices to show that, subject to (6.8), \Q\ > 2. 

By (6.8), there exists a subsequence (m^ : k > 1) along which {CeCrf^^) converges 
to either a or —a. Assume the hrst; the proof is essentially the same in the second 
case. For simplicity of notation, we shall assume that 

lim {aeafj = a. 

m^oo 

By the invariance of under sign change of the conhguration, 

lim /ioo(o-e = 1 I CT/^ = 1 ) = 5(1 -h a), 

m^oo 

lim /ioo(cre = l\af^ = - 1 ) = ^(1 - a). 

m^oo 

Find M such that 

hoo(t^e = 1 = -1) < i(l - la) 

< |(1 + \a) < /ioo(o-e = 1 I CT/^ = 1), m > M. 

We may hnd an increasing subsequence (r^ : m > M) such that 

= 1 = -1) < i(l - |«) 

< 1(1 + |a) < ^lr^{(Te = 1 I = 1), m > M. 

Let /i’*' (respectively, /r“) be a subsequential limit of /ir^(- | = 1) (respectively, 

hr-^(- I = - 1 )), so that 

/i+(cre = 1) > /i"(cre = 1). 

In particular, fi~^ fi~. Since \e — fm\ —)■ cxd as m —)■ 00, the measures /i^ satisfy the 
DLR condition, and therefore they lie in □ 


6.3. Free energy. A boundary condition Bn is a conhguration on the half-edges 
Hn of the planar graph = {Vn, En, Hn), in the notation of Section 6.1. Let 
Zn{a,b,c, Bn) be the partition function of the 1-2 model on with parameters a, 
b, c and boundary condition Bn (as in (4.1), say). The free energy, for given a, b, c 
and boundary conditions {Bn : n > 1), is dehned to be 


(6.9) 


X(a, b, c, (Bn)) 


lim 

n^oo 



log Zn{a,b, C,Bn), 


whenever the limit exists. 
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Proposition 6.4. Let {a,b,c) ^ (0,0,0). The free energy of (6.9) exists and is 
independent of the choice of boundary conditions {Bn) ■ Moreover, up to a smooth 
additive constant, it satisfies 

(6.10) J^(a, b,c) = ^ [[ log P(e'^ e'^) dO dcf, 

47! JJ[0,27r]2 

where P is given in (5.2). 

Proof. The correspondence between 1-2 model configurations on ]HI(j (with the mixed 
boundary condition) and dimer configurations on ]HI„ a was explained in Section 

6.1. It follows that the free energy of that 1-2 model is the same as that of the 
corresponding dimer model. The expression (6.10) follows for that case from a general 
argument used to compute the free energy of this dimer model, given that either the 
spectral curve does not intersect the unit torus, or the intersection is a unique real 
point of multiplicity 2. See Proposition 5.1 and also [23, Thm 3.5] and [3, Thm 1]. 

Next we prove that the free energy of (6.9) is independent of the choice of {Bn). To 
this end, we consider the boxes and illustrated in Figure 5.2. We claim that, 
for any boundary condition on EI)^ (that is, any present/absent configuration on Hn) 
and any 1-2 model configuration on ]HI /_2 (so that the edge-states on P „_2 U iPn -2 are 
given), there exists a configuration on \ (Pn -2 U Pn- 2 ) such that the composite 
configuration is a 1-2 configuration on M/. 

This claim is shown as follows. Consider a given boundary condition on Hn and a 
1-2 configuration on IHI(^_ 2 - The vertex-set Vn\Vn -2 forms a cycle Cn with even length, 
together with two further vertices Wi, W 2 at the left and right corners, see Figure 

5.2. From Cn we select a perfect matching. By considering the various possibilities, 
we may see that it is always possible to allocate states to the two edges between Cn 
and the Wi in such a way that, in the resulting composite configuration, each Wi has 
degree either 1 or 2. 

Let Bn be the free boundary condition, under which no half-edge is present. We 
have that 

\\ogZn{a,b,c,Bn)-\ogZn- 2 {a,b,c,B^n- 2 ) \ < |K\i4-2|^, 

for some K = K{a,b,c) > 0 and all Bn- Divide by |14| and let u —)■ cx) to obtain the 
claim. The theorem follows on noting that the number of boundary configurations 
is 2l'^"'l, and |if„|/|W| —)■ 0 as n —)■ cxd. □ 

7. Proofs of Theorem 3.1(a) and of the limit in Theorem 3.1(c) 

The basic structure of the proof is as follows. As in Section 5, the 1-2 model may 
be represented as a dimer model on a certain decorated graph ]HI„ a derived from 
]HI„. Subject to condition (3.1), the two-edge correlation {o'e<Jf)n of the 1-2 model 
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may be represented in terms of certain cylinder probabilities of the dimer model. 
Using the theory of dimers, these probabilities may be expressed in terms of ratios 
of Pfaffians of block Toeplitz matrices, and a similar representation follows for the 
inhnite-volume two-edge correlation (cxecrj). By Widom’s theorem [37, 38], the limit 
A(a, b, c) := hm|e_/|^oo(c’‘eC’'/)^ exists, and furthermore A is analytic except when the 
spectral curve intersects the unit torus. This identihes the phases of the 1-2 model, 
and they may be identihed as sub/supercritical via the extreme values of (3.3). 

We shall refer to the following condition on the edges e, / of H: 

e and / are midpoints of two NW edges such that 
(7.1) there exists a path tt = 7r(e, /) in AM.n from e to f 

using only horizontal and NW half-edges. 

See Figure 7.1. The principal step in the proof is the following. 



Figure 7.1. A path vr comprising horizontal and NW mid-edges, con¬ 
necting the midpoints of two NW edges e and /. 

Theorem 7.1. Let e,f be two edges satisfying (7.1), and let a>b >0. The limit 
A{a,b,c) := lim|e-/|->.oo(c’'eC’'/)^ exists and is complex analytic in c > 0 except when 
\/c = ^/a — \/b and a/c = y/a -|- \/b. 

Before giving the proof of Theorem 7.1, we explain how to deduce some of the 
claims of Theorem 3.1 (a, c); the exponential rate of convergence in part (c) is proved 
in Section 8. 

Proof of Theorem 3.1(a, c), without the exponential rate. Part (a) holds by Theorem 

6 . 2 . 

Let a > 6 > 0, and let e, / satisfy (7.1). By Theorem 7.1, the function A 3 (-) : = 
A(a, b, ■) is complex analytic on each of the intervals 

6-1= [0,{V^-Vbf), 

C2 = ((a/o - \/&)^, (a/o + \/&)^) , 

C 3 = ((a/o -F Vly, 00 ). 
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(That is to say, for c & Ci considered as a line in the complex plane, A 3 is analytic 
on some open neighbourhood of c.) 

Remark 7.2. Note in passing that, by Remark 4-2, we have A 3 (c) = 0 for c lying 
in the interval S := (a/o^ — 6 ^, + IR). Since S' C C 2 and A 3 is analytic on C 2 , 

we have as implied in part (b) that {ae,<Jf) —)■ 0 as |e — /| —)■ cxd. (This is trivial if 
a = b.) 

We turn to part (c). Consider hrst the interval Ci, and assume a > b. Since 
non-trivial analytic functions have only isolated zeros, it follows that: either A 3 = 0 
on Cl, or A 3 is non-zero except possibly on a set of isolated points of Ci. By (3.3), 
A 3 ( 0 ) = 1 , whence the latter holds. 

By (3.3), (ceCTf) = 1 when a = b = 0 and c = 1. Since A 3 is analytic (and hence 
continuous) on C 3 , there exists a > 0 such that A (a, 6 , c) > a in a small (real) 
neighbourhood of (0, 0,1). Since A depends only on the ratios a : b : c (cf. (2.5)), we 
deduce that, for hxed a, 6 > 0 and sufficiently large c, we have A 3 (c) > a > 0. By 
Theorem 7.1, A 3 is analytic on C 3 , and the claim holds as above. □ 

The remainder of the section is devoted to the proof of Theorem 7.1. We shall 
develop the notation and arguments of Section 6.1. Let be the 1-2 measure on 
]HI(j with the mixed boundary condition of Section 6.1, and let /i™"" := hmn-s .00 h™ 
as after Theorem 6.1. By Theorem 6.2, the 1-2 measure on EI^ satishes —>■ 

Too = as n —)■ cx). 

Let e, / be edges of the hexagonal lattice El satisfying (7.1). Let the path tt of 

(7.1) traverse a total of 2fc — 1 edges and two half-edges, so that vr passes 2k bisector 
edges of the inhnite decorated graph El a. We denote this set of bisector edges by 

(7.2) B={bi = {ui,Vi):t = l,2,...,2k}, 
where Vi & tt. 

Our target is to represent {(TeCTf) as the Pfaffian of a truncated block Toeplitz 
matrix, as inspired by [20, Sect. 4.7]. A principal difference between [20] and the 
current work is that, whereas bipartite graphs are considered there and the determi¬ 
nants of weighted adjacency matrices are computed, in the current setting the graph 
is non-bipartite and we will compute Pfaffians. 

To EIa we assign a clockwise odd orientation as in Figure 5.3: the hgure shows a 
clockwise odd orientation of Eli a, embedded in a 1 x 1 torus, that lifts to a clockwise 
odd orientation of EIa. As in (5.1), a horizontal (respectively, NW, NE) bisector edge 
of EIa is assigned weight a (respectively, b, c), and all the other edges are assigned 
weight 1. The bisector edges Pi = {ui,Vi) are oriented in such a way that each pi is 
oriented from Ui to Vi in this clockwise-odd orientation. 

Let Kn be the Kasteleyn matrix of EI„ a (as in Section 5.2 and [25]), and let |t| 
denote the index of the row and column of corresponding to the vertex v. Assume 
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that luji = |Mj| + 1 for 1 < i < 2/c, and furthermore that 

(7.3) |ki| < luil < \U 2 \ < 1 ^ 2 ! < ■ ■ ■ < \u 2 k\ < \v 2 k\- 

Let K~^ be the inhnite matrix whose entries are the limits of the entries of K~^ as 
n —)■ 00 . The existence of K~^ may be proved by an explicit diagonalization of 
using periodicity, as in [9, Sect. 7] and [24]. 

We now construct the modihed Kasteleyn matrix Ki{z,w) of Hi^a by multiplying 
the corresponding entries in its Kasteleyn matrix by z or z~^ (respectively, w or 
according to the manner in which the edge crosses one of the two homology 
cycles 7 j,, '^y indicated in Figure 5.3. As remarked in Section 5.2, the characteristic 
polynomial P{z, w) = det Ki{z, w) is the function /(a, 6, c; tc, z) of (5.2), see also [25, 
Lemma 9]. The intersection of the spectral curve P{z,w) = 0 and the unit torus 
is given by Proposition 5.1. 

Consider the toroidal graph EI^ a. Let 'yy be homology cycles of the torus, 
which for dehniteness we take to be shortest cycles composed of unions of boundary 
segments of fundamental domains as in Figure 5.3. Let Kn{z,w) be the modihed 
Kasteleyn matrix of IHI^^a. 

For I C {1 ,2,..., 2k}, let Mj be the event that every bi with i G / is present in 
the dimer conhguration. Assume n is sufficiently large that 

(7.4) for 1 <i < 2k, the edge intersects neither nor y^,. 

For n even, as in (6.4), 

(7.5) 

71oo(A 7/) = lim 7 i„(M,) 

n^oo 

^ J. -Pf jFny(i, 1) + Pf -1) + Pf k^A-i, 1) + Pf iL^,j(-i, -1) 

^ -PfiF,,(l,l) + PfiF„(l,-l) + PfiF,(-l,l) + PfiF,,(-l,-l) ’ 

where 

(7.6) j,= 

iei 

and Aj denotes the submatrix of the matrix A after deletion of rows and columns 
corresponding to {ui,Vi : i E 1} (see [25, Thm 0.1]). Note that /4oo(Af0) = 1. 

The limit of (7.5) can be viewed as follows. Each monomial in the expansion of 
Pf Kn{z,w), z,w E { — 1,1}, corresponds to the product of edge-weights of a dimer 
conhguration, with possibly negative sign. The coefficients in the linear combination 
of Pf Kn{l, 1), Pf Kn{l, —1), Pf Kn{—1, 1), aud Pf Kn{—1, —1) are chosen in such a 
way that the products of edge-weights of diherent dimer conhgurations correspond 
to monomials of the same sign. The numerator of (7.5) is the sum over dimer 
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configurations containing every bi, i G I; this can be computed by the corresponding 
sum of monomials in the expansion of the denominator. Under (7.4), : = 

[Kn{z,w)]ui^vi is independent of z,w E { — 1,1}. Since each bi is oriented from m, to 
Vi, we have that [Kn]ui,vi = c, whence Jj = 


Lemma 7.3. Let A be a 2m x 2m invertible, anti-symmetric matrix, and let L C 
{1,2,..., 2m} he a nonempty even subset. Let he the submatrix of A obtained by 
deleting the rows and columns indexed by elements in L, and let Af^ be the submatrix 
of A~^ with rows and columns indexed by elements in L. Then 

[-ifiApfAL = Pf (7l)Pf where S{L) = 

l&L 


Proof. See, for example, [5, Lemma A.2]. 


□ 


The conclusion of Lemma 7.3 holds also when L = 0, subject to the convention 
that Pf {Aff^) = 1. 

Returning to (7.5), take L = {|Mj|, |ni| : i G /}, so that (—1)'^^'^^ = (—1)1^1, by the 
choices before (7.3). When the spectral curve does not intersect the unit torus T^, 
by Lemma 7.3, 


lim 

n^oo 


PfKnj{z,w) 

PfKn{z,w) 


lim (—l)l^lpf K.^]{z,w) 

n^oo ’ 

{-lf\pfKf\ 


where the limit is independent of z,w E ( — 1,1}; see [24, Lemma 4.8] for a proof of 
the existence of the limits of the entries of K~^. By (7.5)-(7.6), 

(7.7) fi^{Mi) = {-cf^PfKY\ 


We shall make use of the following elementary lemma, the proof of which is omitted. 


Lemma 7.4. Let S be a random subset of the finite nonempty set B. The probability 
generating function (pgf) G{x) = E(xl'^l) satisfies 

G(l + A) = ^Al^lp(AD J), AeM. 

ICB 

Let B be the set of bisector edges along vr (see (7.2)), and let S be the subset of 
such edges that are present in the dimer conhguration. By (6.6), 

(aeCT/) = G(-l), 

where G is the pgf of 151 under the measure poo- By (7.7) and Lemma 7.4, 

(7.8) {a,a,) = ^(-2)l'Lco(A7;) = ^(2c)l'lpf AT'. 

ICB ICB 
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This may be recognized as the Pfaffian of a certain matrix defined as follows. 
Let Li(A) be the 2x2 matrix 


l'.(A) = 1 ) . 

and let Y 2 k{)^) be the 4k x 4k block diagonal matrix with diagonal 2x2 blocks equal to 
yi(A). More precisely, T 2 fc(A) has rows and columns indexed Ui, Vi, U 2 , V 2 , ■ ■ ■, U 2 k,V 2 k, 
and 



/A (A) 

0 

■ ° \ 

Y2k{\) = 

0 

A(A) ■■ 

0 


V 0 

0 

■ A(A)/ 


Lemma 7.5. We have that 


(7.9) {aeaf) = Pi[Y2kil) + 2cKy^\ 
where 14 = {ui,Vi,U 2 ,V 2 ,.. .,U 2 k,V 2 k}- 

Proof. It suffices by (7.8) that 

(7.10) Pf[Y2k{l) + A] = J2PiAi, 

ICB 


where A = (ajj) is a x 4k anti-symmetric matrix with consecutive pairs of 
rows/columns indexed by the set B = {1,2,..., 2fc}, and Ai is the submatrix of 
A with pairs of rows and columns indexed by / C S. 

Let G = {V,E) be the complete graph with vertex-set V = {1,2,..., 4fc}, and 
recall that 


(7.11) PiA = 5^sgn(7r^) J] a,j, 

/isn 

i<j 

(see [18, 35]), where fl is the set of perfect matchings of G, and the permutation 
vr^ e 54^ is given by 


A 2 3 4 ■■■ 4fc-l 4k\ 

V* jl 4 4 ■ ■ ■ i2k j2k) 


where fi = {(4,4) : 1 < r < 2k}, 4 < 4 < • • • < and 4 < 4- 
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By (7.11), 

(7.13) Pf \y2kiX) + A] = sgn(7r^) n 

Aisn (*j)eAt 

i<j 

= 5^fesgn(7r^) Yl [^ 2 fc(l)]i,i n 

Kcv Aten (Li)eM 

ieK, i<j i^K, i<j 

The penultimate product is 0 unless every i G 77 is odd and satisfies {i,i + 1) G /i. 
Therefore, with J = + 1 ) C i?, 

Pf [Y2k{l) + A] = Y^ sgn(7r^) JJ atj 

JCB /jSII i<j 

ij^(2J-l)U(2J) 

JCB 

as required for (7.10). □ 

Recall that a matrix is Toeplitz if every descending diagonal is constant, and a 
block matrix is block Toeplitz if each block is Toeplitz and every descending diagonal 
of blocks is constant. Now, Y 2 k (1) + 2cKy^ is a truncated block Toeplitz matrix each 
block of which has size 4x4. We propose to use Widom’s formula (see Theorem 7.8) 
to study the limit of its determinant as 7 —?• cx). In this limit, the matrix becomes 
an infinite block Toeplitz matrix T('^) with symbol given by 

(j){z, e*®) dO 

where 

(7.15) 0(^,e*^) = F2(1) + 2ci7r'(^,e*')(n4), 

and ^( 1 - 4 ) denotes the 4x4 submatrix of the matrix A with rows and columns indexed 
by Mi,ni,M 2 ,T 2 as in Figure 5.3. This follows by the explicit calculation 

-I p27t p27r 

(7.16) [K-X,, = —j j u,veV^, 

where v' is the translation of v to the same fundamental domain as u, and k is 
the number of fundamental domains traversed in moving from v' to v, with sign 
depending on the direction of the move. When 7 7 ^ 0, (7.16) is the kth Fourier 
coefficient of the symbol (7.14). See [24, Sect. 4] for a similar computation. 


(7.14) 


r2w 


V'(^) = ^ 
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Lemma 7.6. Let z G C with \z\ = 1. When the spectral curve does not intersect the 
unit torus we have that diei'ip{z) = 1. 

Proof. Let ]HIm,n,A be the toroidal graph comprising m x n fundamental domains (a 
fundamental domain is drawn in Figure 5.3). We can think of ]HIm,n,A as the quotient 
graph of Ha under the action of mZ x nZ. To Hm,n,A we allocate the clockwise-odd 
orientation of Figure 5.3. 

Let Km,n{z, w) be the corresponding modihed Kasteleyn matrix. Assume the cycle 
73 ; (respectively, 7 ^) crosses 2n (respectively, 2 m) edges, whose weights are multiplied 
by X or z~^ (respectively, w or w~^), depending on their orientations. Note that 
= K,. 

The toroidal graph Hi a is a line of n copies of the graph of Figure 5.3, aligned 
parallel to 73 ,. It contains 2n (bisector) edges with weight c, of which we select 
two, denoted ei, 62 , lying in the same fundamental domain. Let Hi^^ be the 
oriented graph obtained from Hi^^i^a by reversing the orientations of ei and 62 , and 
let KIj^{z,w) be the modified Kasteleyn matrix of Hi„^. 

Let X(A) be the 4n x 4n matrix 

/Ki(A) 0 0\ 

A(A) = 0 Ki(A) 0 . 

\ 0 0 0 / 


Since ei and 62 have weight c, 

^ detKl^{z,w) _ det[Kl^{z, w) - Ki^njz, w) + Ki^njz, w)] 

det Ki^niz,w) detKi^niz,w) 

_ det [X (-2c) -h Ki^n{z, w)] 
det Ki^n{z, w) 

= det [X{-2c)Kf^{z,w) + I] 

= det [K2(-2 c)A:/^(x, w)(i:4) + /] 

= det [2cKf^{z,w)^i,4) + I2 ( 1 )] , 


for w = ±1, since Y2{1)Y2{—1) = I and det > 2 ( 1 ) = 1- Here, iFi/(x, w)(i: 4 ) is the 
submatrix of Kf^{z, w) comprising the rows and columns indexed by the four vertices 
incident with the Cj, see Figure 5.3. 

By an explicit diagonalization of as in [9, Sect. 7], for any two vertices u, v 
in the same fundamental domain, the limit 


1 



d9 


(7.18) 


bm [K^ f{z,w)]u,v 

n^oo 
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exists and is independent of the choice of w = ±1. The proof of the next lemma is 
deferred nntil the cnrrent proof is completed. 


Lemma 7.7. For z E C, the limit 


(7.19) 6{z,w) 

satisfies 6{z, —1) = 6{z, 1) = ±1. 


lim 

n^oo 


det Kl^{z,w) 
detKi^n(,z,w) 


We dednce that 


det ^p{z) = lim det [> 2 ( 1 ) + 2cK^^{z,l)n.4)\ by (7.14), (7.15), (7.18) 

= 6{z, 1) = ±1 by (7.17) and Lemma 7.7. 

Setting z = 1, we have by (7.19) that det "0(1) > 0 since it is the limit of a ratio of 
determinants of two anti-symmetric matrices. By (7.14) and the forthcoming (7.25), 
ip is continnous on the nnit circle when the spectral cnrve does not intersect T^, and 
the claim follows. (Note that (7.25) is a general fact whose proof does not depend 
on Lemma 7.6.) Therefore, '4){z) = 1 for \z\ = 1. □ 


Proof of Lemma 7.7. By (7.17)-(7.18), 5 ( 2 ;,—1) = S{z,l) =: S{z), say. We claim 
that 


(7.20) det Kf^{z, —1) = det Ki^n{z, 1), det 1) = det Ki^n{z, —1). 

By (7.19)-(7.20), 6{z, —1) = l/5(^, 1), so that 6{z) = ±1 as claimed. 

We prove (7.20) next. Each non-vanishing term in the expansion of det Kl^{z, —1) 
and detiLi^„( 2 :, 1 ) corresponds to a cycle conhgnration on that is, a conhgn- 

ration of cycles and donbled edges in which each vertex has two incident edges. Let 
C be an oriented cycle of ^ ^ viewed as an nnoriented graph, and recall that Ci 
and 62 have opposite orientations in Hi „ a and H^ It snffices that C contribntes 
the same sign on both sides of the left eqnation of (7.20). Let c(C) (respectively, 
w{C)) be the nnmber of c-type (respectively, tc-type) edges crossed by C. By a 
consideration of parity, c{C) is even if and only if w{C) is even, and in this case 
C contribntes the same sign. We claim that w{C) = 1 if c(C) = 1. It is standard 
that C is either contractible or essential, and that C, if essential, has homology type 
±1 in the direction 7 ^,. Therefore, w{C) = 1, and the claim follows. The second 
eqnation of (7.20) follows similarly. □ 


We remind the reader of Widom’s theorem. 
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Theorem 7.8 (Widom [37, 38]). Let Tm(0 be a finite block Toeplitz matrix with 
given symbol ^ and m x m blocks. Assume 


where 


E ii&ii + E i*=i ■ II& 


< oc, 


k=—oo 
J0\ 


\k=—oo 


det ^^ 0, 7 ^ Ao<e< 2 ,r arg det ^(e*®) = 0, 

ZTT 

denotes Hilbert-Schmidt norm, is the kth Fourier coefficient of and 


^A„<,<,,argdet5(e") = ^ 


f <iC J:logdet5(C). 

'|CI=1 ^8 


Then 


where 

(7.21) 

(7.22) 


lim 

m^oo 


detTmiO 




G(0=exp|^^ logdet^(e*®)d0| , 
E{0 = det [T(OT(r')] , 


where Ti^C) is the semi-infinite Toeplitz matrix with symbol and the last det refers 
to the determinant defined for operators on Hilbert space differing from the identity 
by an operator of trace class. 


We note that, when the spectral curve does not intersect the unit torus, if given 
by (7.14) is a smooth matrix-valued function on the unit circle, whence 

ll^fcll < Ca\ 

for some C > 0, a G (0,1), and all fc > 0, where ifk is the /cth Fourier coefficient of 
if. 

Proof of Theorem 7.1. This holds as in the proofs of [24, Lemmas 4.4-4.7], and the 
full details are omitted. Here is an outline. 

The symbol i/i is a 4 x 4 matrix-valued function. Let n be the number of NW edges 
in the path of (3.1) connecting e and /. By the computations of Section 7 (see (7.9) 
and (7.14)), 

(7.23) {(Tecr/)^ = det T„(V>), 

where T^fif) is a truncated block Toeplitz matrix consisting of the hrst nxn blocks 
of an inhnite block Toeplitz matrix T{if) with symbol if (so that T^fif) is a 4n x 4n 
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matrix). By Lemma 7.6, when the spectral curve does not intersect the unit torus 
T^, we have = 1, where G{'ijj) is given in (7.21). By Theorem 7.8, 

(7.24) A(a,b,c) = lim detTn(V’) = 

n^oo 


where is given in (7.22). 

Non-analyticity of A may arise only as follows. One may write 


(7.25) 


[ATi ^{z,w)]ij 


Qij{z,w) 

P{z,w) 


where Qij{z, w) is a Laurent polynomial in z, w derived in terms of certain cofactors 
of Ki{z, w), and P{z, w) = det Ki{z, w) is the characteristic polynomial of the dimer 
model on It follows that A is analytic when P{z,w) has no zeros on the unit 

torus T^. The last occurs only under the condition of Proposition 5.1, and the claim 
follows. □ 


8. Proof of exponential convergence in Theorem 3.1(c) 

We develop the method of proof of Widom’s formula. Theorem 7.8, see [37, 38]. 
For a positive integer r, let Ar H be the Banach algebra of r x r matrix-valued 
functions on the unit circle under the norm 

oo / oo 

II0II = 1^1 ■ 

k=—oo \k=—oo 

where 0^ is the fcth Fourier coefficient of the matrix-valued function 0. Note that 
the trigonometric polynomials are dense in Ar H Kr- 
The main theorem of this section is as follows. 

Theorem 8.1. Let T{'ijj) be a semi-infinite block Toeplitz matrix with symbol 
where is an r x r matrix-valued, -function on the unit circle. Let Tn{'f) be the 
truncated block Toeplitz matrix consisting of the first n x n blocks of In the 

limit as n ^ oo, detT„('0) converges to its limit exponentially fast. 

We recall (7.23) from the last section, and note that "0 is a matrix-valued function 
in Ai n 7^4. Let be the Hankel matrix with symbol fj, 

Ei'f) — ('0*+t+l)o<ij<oo’ 

and write 

( 8 . 1 ) ^{z) =^p{z-^). 

Lemma 8.2. Let e, f be NW edges of the hexagonal lattice satisfying (3.1). In the 
representation (7.23) of the edge-edge correlation fje,Of)n, the symbol xf is C°° on 
the unit circle whenever the spectral curve does not intersect the unit torus. 
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Proof. We recall from Proposition 5.1 that, when a > 6, c > 0 and 

^/a Vb + \/c, a/c 7 ^ y/a + Vb, 


the characteristic polynomial P{z,w) has no zeros on the unit torus. As in Section 
7, '0(0 is a matrix-valued function dehned on the unit circle, each entry of which 
has the form 


1 


f27r 


0(000 


d(j), 


P(0e*0 

where P{z,w) is the characteristic polynomial, and Q{z,w) is a Laurent polynomial. 
When P{z, w) has no zeros on the unit torus, "0(0 is a function on the unit circle, 
and the nth Fourier coefficient of "0(0 decays exponentially to 0 as |n| —)■ cx). □ 


Assume first that the operator T(0) is invertible as an operator on I 2 
seqnences of r-vectors. This is equivalent to assuming that the matrix-valued 
function 0 has a factorization of the form 

( 8 . 2 ) 0 = 0 + 0 _, 

where the 0± are invertible in fl Kr, and the 0^^ (respectively, 0^^) have Fourier 
coefficients that vanish for negative (respectively, positive) indices. As in [38, Sect. 
3], we have 

(8.3) ‘ awnt ^ = '‘«‘(^ - PnH ( m ( P -) PnT { r +) Pn ), 

where PnA is the submatrix of A consisting of its first nr rows, and BP^ is the 
submatrix of B consisting of its hrst nr columns. Recall that G is given by (7.21). 
Now we dehne operator norms, and discuss inequalities regarding these norms. 

Definition 8.3. For a compact operator A on a Hilbert space, and 1 < p < 00 , 
let ||A||p denote the p-norm of the eigenvalue-sequence of {A*A) 2 , where A* is the 
conjugate transpose of A. The oo-norm is the usual operator norm and is so defined 
even if A is not eompact. The 2-norm is the Hilbert-Schmidt norm, and the 1-norm 
is the trace-norm. The set of compact operators with finite p-norm is denoted by 
T>p. Thus <Pi is the set of operators of trace class; T >2 is the set of Hilbert-Sehmidt 
operators; and <Poo is the set of compact operators. 

As in [38, Sect. 2], we have the following lemma 

Lemma 8.4. 

(a) Let 1 < p < 00 . If A ^ <Pp and B,C & <Poc, then BAG G <Pp. Moreover, 

\\BAG\\p<\\A\\p\\B\U\C\U 
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(b) If A, B G <p 2 , then AB G Moreover, 

\\AB\\, < \\AUBh. 

Lemma 8.5. Assume thatif{() is a C°° rxr matrix-valued function on the unit circle 
with exponential decaying Fourier coefficients. Assume also that has a factorization 
given by (8.2). Then all C°° rxr matrix-valued functions 

on the unit circle with exponential decaying Fourier coefficients. 

Proof. By the arguments of [38, p. 10], when T('^) is invertible, is a matrix-valued 

function whose sequence of non-negative Fourier coefficients is T“^('^)(/, 0, 0,...). 
By [33, Thm 1.3], the entries in the hrst column of T{gp)~^ decay exponentially as 
the entry moves away from the diagonal, whence the nth Fourier coefficient of 
decays to zero exponentially as |n| —)■ oo. 

The exponential decay of Fourier coefficients and smoothness of '0+^; 
follow. □ 

By explicit computations (see [38]), 

Moreover, 

1 

< \\Hm\2 ■ l|ff(F‘)(a - /)||2 ■ i|T(^;')iico. 

Lemma 8.6. Assume that ip is a rxr matrix-valued, C°° function on the unit circle 
with exponential decaying Fourier coefficients, and assume that has a factorization 
given by (8.2). Let Tr(y4) be the trace of A. Then there exists 0 < /3i < 1 such that 

(8.4) ||T(i^;')|U<oo, 

(8.5) ||W )||2 = [^TiH*{^p)H{^P))]-^ < oo, 

(8.6) \\H{^Z^){Pn - I)\\2 = TT[{Pn " I)H*{^Z^)H{^Z^){P^ - /)] ^ < /Si". 

Proof. Note that ||T('0^^)||oo is exactly the maximal singular value of TpipZ)^). Note 
also that T('^+)“^ is an upper triangular matrix whose off-diagonal entries decay 
exponentially fast. Let 
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It is not hard to check that the off-diagonal entries of A decay exponentially fast 
with respect to their distance to the diagonal. Moreover, is the maximal 

eigenvalue oi A. It is well known that 

OO 

||T(V';^)||L < max^ |aij| < oo, 

^ i=i 

and (8.4) follows. The inequality (8.5) is obtained since the diagonal entries of 
decay exponentially when moving down the diagonal. Equation (8.6) is 
similar. □ 


Now, 

(8.7) 

< \\PnH{^)H{^-^){Pn - J)ll^ + 1|(P„ - 

< ||P„|U ■ ||P(l^)P(^-')(Pn - /)||l + ||(Pn - I)Hm\2 ■ II W"')l|2 

< ||P(l^)||2 • II W-')(Pn - /)||2 + \\{Pn " /)W)||2 ' || W-')||2 
</^ 2 , 


where 0 < /92 < 1- The last inequality holds because the {i,j) entries of P('^) and 
decay exponentially in i + j. 

The following cases may occur: 

I. J — PI{'ip)PI{'ip~^) is invertible, 

II. I — is not invertible. 

Assume Case I occurs. By the formula of [13, p. 116], 


( 8 . 8 ) 


det(/ - PnH(^)H(^z')PnTH’l')Pn) - det(/ - 

< g||R(y)ri(y-blli(eQn _ 


where 

Qn ■ 




and 0 < < 1. 

Assume Case II occurs. We follow the approach of [13, pp. 116-117]. When I — 
is not invertible, the point 1 is an isolated eigenvalue of hnite type for 
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H{'ip)H{'ip ^). Let P be the corresponding Riesz projection, and put Hi 
H 2 = KerP, so that Hi is hnite-dimensionah For simplicity, we write 


Fn := PnH{^/J)H{^/JZ^)P^T{^/J-^)Pn, 

A := H(^)H(^-^). 


With respect to the decomposition H = Hi (B H 2 , we have 


P„ = 


-'^12 
Rp) Ri 
-'^21 -'^22 


and it follows as in (8.7) that 


A 


f All 0 \ 

\ 0 A 22 J ’ 


\\Fn-A\\2<\\Fn-A\\i<P2, 


ImP, 


for some 0 < /?4 < 1. Moreover, 

{Fn-Ay{Fr,-A) = 


Li 


where 


Li = _ Aiiy{Ky> - Rn) + {KyirKyi\ 


Zn) 


iA\* t^A) 


L2 = {Kt^yKy^> + {K)^> - A 22 y{Kyy - 7I22) 


An) 


(n) 


-(n) 

‘■22 


A 


11 


+ Tr 


(n)^ 
21 y 


K 


Therefore, 

||F„ - Alii = ((/<■£> - 

+ IV ((A'£> - A^^riAt’ - mZ) + Ti' {(A 
= II A'!)' - AnWl + ll/h’di) + II a!;>||) + ||/4> - A: 


(n) 

21 


(n)^ 
12 ; 




(n) 

12 


22 || 2 - 


Hence, 

iia 1;> - Aiiih < ft”, 

IIA'g'lh < ft”, 

IIAt'lh < ft”. 

As in [13, pp. 116-117], with each Ij an identity matrix of suitable size, 
det(/ - P„) = det (a - Kif + kA\I 2 - det (/2 - H^), 

and 

lim det(J 2 — ^ 22 ^) = det(/ — ^ 22 )- 

-T)_ ^ 


36 


GEOFFREY R. GRIMMETT AND ZHONGYANG LI 


Moreover, 


h - 


("■)/ 




ih - ^ll) 


< -An\\i + \\K- 


W| 
12 I 


\h-K. 


Wi 


22 l|co 


\K. 


(»^)i 


21 


for some 0 < /Sg < 1- Since Ii — + K^{l 2 — K^) is a finite matrix, all 

the norms are eqnivalent, we have 

det (Ji - + K^;^{I 2 - - det(/i - Aln)| < 

for some 0 < /^y < 1. 

Since (by Lemma 7.6) G{'ip) = 1, we obtain that, when Ti^ip) is invertible, (cJeCT/-)^ 
converges to its limit exponentially fast in the limit as \e — f \ —)■ oo. This completes 
the consideration of Case II. 


Assume now that T{'ijj) is not invertible. Since T{'ijj) is Fredholm of index 0, by 
[38], there exists (j) with only finitely many non-zero Fourier coefficients, such that 
T{'il> + ecf)) is invertible for sufficiently small e 7 ^ 0. Therefore, 

„'i?i = '*'* + £^)-‘)]. 

for e belonging to some punctured disk with centre 0. Using the same arguments as 
above we obtain that, for sufficiently small |e| > 0 , 


(8.9) 


detTn{i! + 60 ) 
G(0 + 60)^+! 


det[T(0e0)T((0 + 60) ^)] 


</5 


n 

5 


for some 0 < 05 < 1 independent of 6 . This is obtained by first expressing the 
ratio on the left side of (8.9) as in (8.3), then follow the previous argument. Note 
that 05 depends on the exponential decay rate of (0 -f 60 )^ as fc —)■ cxd, which can 
be made universal for e with | 6 | sufficiently small. Since the component functions 
on the left side of (8.9) are analytic in e on some disk {6 G C : |6| < cq}, by the 
maximum principle, (8.9) holds for 6 = 0 also. Hence ((Jea/)^ converges to its limit 
exponentially fast. The proof of Theorem 8.1 is complete. 


9. Alternative proof of Theorem 3.1(b) 

We find it slightly more convenient to work here with horizontal edges rather than 
NW edges, and there is no essential difference in the proof. Let e, / G E be horizontal 
edges. There exists a horizontal edge g E K such that: e and g (respectively, / and 
g) are connected by a path l^g (respectively, Ifg) of AH comprising only horizontal 
and NW (respectively, NE) edges, and the unique common edge of leg and Ifg is g. 
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Let m + 1 (respectively, n + 1) be the number of horizontal edges in leg (respectively, 
Ifg). Write m An = min{m, n} and mV n = max{m, n}. 


Theorem 9.1. Assume the parameters of the 1-2 model satisfy 


(9.1) a >b > 0, 

(a) The limit 

exists, and 


C > 0, y/a 7 ^ Vb -t- y/c, y/c ^ yfo, + Vb. 
A{a,b,c) := lim (aeCX/)^ 

m,n^oo 


(9.2) \{aeaf)^-A\<Ca^^^, 

where C > 0 , a G ( 0 , 1 ) are constants independent of m, n. 

(b) If, in addition, y/a < y/b + y/c and y/c < y/a + y/b, then C > 0, a G (0,1) 
may he chosen such that 

(9.3) |{crecr/)| < 


Proof, (a) By a computation similar to that leading to (7.9), {aeUf/ may be ex¬ 
pressed in the form 


{aeay/ 



- 

T„(V>2) 


1 


where -0* is a 4 x 4 matrix-valued function (see (7.14)), T/f)/) is a 4r x 4r truncated 
block Toeplitz matrix with symbol "0,, and Am,n (respectively, Bn,m) is a 4m x 4n 
(respectively, 4n x 4m) matrix with entries satisfying 

\A,^Ahj)\ < C/A. \Bn,AAj)\ < C/A. 


where C > 0, 01 G (0,1) are constants independent of i, j. 

Let 01 , 02 be two matrix-valued functions dehned on the unit circle with only 
hnitely many non-vanishing Fourier coefficients, such that T(0i-f-e0i) and T( 02 +e 02 ) 
are invertible for complex e 7 ^ 0 with sufficiently small modulus (recall (8.1)). Let 
e > 0 be given accordingly, and let 


Ae = Ti + e0i, i = 1, 2. 

Since T{/jA is invertible, we have (as in Section 8 ) that 

where the 0i,e,± are invertible in A^PKi, and //}+ (respectively, /fl-) have Fourier 
coefficients that vanish for negative (respectively, positive) indices. 

Subject to (9.1), by Proposition 5.1 the spectral curve has no zeros on the unit 
torus. As explained after Theorem 7.8, the entries of the Toeplitz matrix (Fourier 
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coefficients of a smooth function on the unit circle) decay exponentially as their 
distances from the diagonal go to inhnity. 

By computations similar to those of [38, p. 8], 


(9,4) 






^m,n 

Bn,m Tji(lp2^e 


5n,mT^(V’l,e,-)^m(V’l,,,+) Tn{^2,e)Tn{i’2,e,-)Tn{i^2,e,+) 


= I 


4m+4n 




m,n? 


where Ir is the r x r identity matrix and 


^ ^ -A^^nTn{ihi-)Tn{iJ2X+) 

PnH {ij2,e) H PnTn{i^2,l-) 


We now take determinants of (9.4). As in [38, p. 8], 


G(V.,.) 


m +1 ’ 


whence 

(9.6) 

Let 


G'(V. 1 ,,)-+ 1 G'(V> 2 ,J "+1 






where A and B are inhnite matrices obtained as the limits of Am,n and Bn,m, as 
m, n —)■ oo. Note that S' is a trace-class operator. Therefore, det(/ — S) is well- 
dehned and complex analytic in a, b, c, whenever the entries of S are analytic in 
a,b,c (see [14], and also [13, Lemma 3.1] and [24, Lemma 4.6]). Then, 


\\S^,n - ^lll < \\Am,nTr.{^P^l_)Tn{^P^l^) - AT{^P^l_)T 

BT{i;Yl_)Tii; 


|5n,mTm(V’lj_)Tm(V'lj,+) 
L -1 


1 

-1 

li^iA 


1 

-1 




+dli 
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Using similar arguments as in (8.7), we can show that 

where /32 G (0,1) is a constant independent of m, n, e. Moreover, 

subject to (9.1). Here, i/’. e A. A e (0,1) 

independent of i,j,e. Therefore, 

j)l<Ci/3U. 

j)l< CiPi*’, 


-2"l,-)(bj) < Cif3l^^, 

|sr(v>r,L)J’(V’r,.‘,- 

\l,+){hj)\< Cil3l^\ 


T.V)(h7)l< 



where Ci > 0 ^ (3^ E (0,1) are constants independent of i, j, e. 
Hence, 


P^AT(ixl^)T(ixl_)P„ - AT(ixlJT(iX 


',-1 






mAn 


where C 2 > 0 , (3^ E (0,1) are constants independent of i,j, e. 
Moreover, 


and 


PmAT{^|J^l_)PnTii^2,l-)Pn " Pm AT i^2,l+) Pn 

= \\Pr^ATi^l^^l_)iPn - I)Ti^lP^l^)P^\\^ 

< \\PmAT{^lP^l_){P^ - /)||i ■ ||T(l72-,V)Pn||oo 

< C,/3A 


P^AP^T{^lj^l_)PnTii^2,l-)Pn - PmAT{^lj^l_)PnTii^2,l+)Pn 
= \\PmA{Pn - I)T{'llj^l_)PnT{'lp2l+)Pn\\l 

< \\PmA{P^ - /)||i ■ ||T(l72-:l,_)PnT(l^2-,V)^n||oo 

< 


where C 3 > 0, (3q E (0,1) are constants independent of m,n,e. 
We consider the following cases: 

I. / — S' is invertible, 

II. I — S is not invertible. 


are 
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If Case I occurs, following arguments similar to those of Section 8 , 

(9.7) |det(J - Sm,n) - det(/ - ^)| < 

where C 4 > 0, (Ij E (0,1) are constants independent of m, n, e. If Case II occurs, by 
considering the Riesz projection and following similar arguments, 

(9.8) |det(J - - det(/ - R)| < 

where C 5 > 0 , /Sg ^ ( 0 , 1 ) are constants independent of m,n,e. 

Letting e —)■ 0 and using analyticity (as in Section 8 ) subject to (9.1), we have 

A := lim (ciecr/)^ = det(/ - S')! 

m,n^oo ' 


by (9.5), (9.7), (9.8) and the fact that Gi'ihi) = 1 (the last follows by Lemma 7.6 and 
(7.21)). Moreover, (10.2) holds. 

(b) The above argument applies also when m is hxed and n —)■ 00 . In this case, 
we replace S in (9.6) by 








where Am,oo, Boo,m are matrices obtained from Am,n, Bn,m by letting n ^ 00 . Now, 
det(/ — Sm) exists since Sm is a trace-class operator, and 

(9.9) lim(creCr/)^ = det(/-S'm)| m > 0, 


as above. We claim that there exists C > 0, a G (0,1), independent of m,n, such 
that 


(9.10) 


(cTeCr/)^ - lim (cTeCr/)" 


< C'a”, m, n > 0 . 


To show (9.10), hrst, following the same computations as above, we obtain 


\\S^,n-S^h<C-jl3;, 

for constants Cy > 0, /I 9 G (0,1) independent of m, n, e. 

We consider the following cases: 

1. J — S' is invertible, 

11. I — S is not invertible. 

In Case I, for sufficiently large m, I — Sm is also invertible. Following the arguments 
of [13, pp. 115-116], we obtain 

|det(/ - Sm,n) - det(/ - Sm) \ < CsPiQ, 

for Cs > 0, (3io G (0,1) independent of m, n, e. 
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In Case II, the point 1 is an isolated eigenvalue of finite type for S. Let P be the 
corresponding Riesz projection, and put Hi = ImP, H 2 = KerP, so that Hi is hnite 
dimensional. With respect to the decomposition, we have 


S' — 


T^im,n) T^{m,ny 
-'^11 -'^12 


5™, = 


T^{m) T^(my 
-'^11 -'^12 

Rif Kf 


We now follow the argument of [13, pp. 115-116] and the proof in Section 8, to obtain 

(9.10) . 

Fix m, and note that is an operator of trace class. By [24, Thm 4.6] (see 

also [14]) and (9.9), hm„^oo(<Te<T/)^ = det(/ — 5'm)|^_Q is analytic in a, b, c subject to 

(9.11) a, 6, c>0, \fa<\fhP\fc, y/h<\fap\fc, y/c < y/a + Vb. 

By Remark 4.2, hm„^oo(<Te(T/)^ = 0 when 

(9.12) a, 6, c>0, < 6^ + c^, b‘^<a^ + c^, c^<a^ + b‘^. 

The set of (a, b, c) satisfying (9.12) is a subset of the (connected) subset of given 
by (9.11), and it follows by analyticity that 

(9.13) lim (cTeCT/)^ = 0, m > 0, 

n—^CXD 


subject to (9.11). 

By the arguments that led to (9.10), there exists C > 0, a G (0,1), independent 
of m,n, such that, subject to (9.11), 


(9.14) 




lim (cTeCr/)^ 

m^oo 


< Ca^, 


m,n > 0, 


and, in addition, 

(9.15) lim = 0) Li > 0, 

m^oo 


We combine (9.10)-(9.15) to obtain 
which implies (10.3) with amended C, a. 

□ 


10. Periodic 1-2 models 

10.1. Two-edge correlation for periodic models. Some of the previous results 
may be extended to certain periodic models, as explained next. Since each edge 
of El touches exactly one white vertex, a 1-2 model may be specihed by assigning 
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a parameter-vector 6^, Cw) to each white vertex w. For k,l E N, we call the 
ensuing model k x I periodic if 

(oiiij biuf Cu,) T~i ^2^1 

where the maps r* are illustrated in Figure 2.3. 

By the techniques of Sections 7, if the parameter-vectors of a periodic model 
are such that the associated spectral curve does not intersect the unit torus, then 
the entries of the corresponding inverse Kasteleyn matrix K~^{v,w) converge to 0 
exponentially as |n — wl —)■ oo. Following the procedure of Section 9, we obtain the 
following (in which the notation of Theorem 9.1 has been adopted). 

Theorem 10.1. Let k,l > 1. Assume the 1-2 model is kxl periodic, and the spectral 
curve does not intersect the unit torus. 

(a) The limit K := lim^ „_,.oo(cre(J/)^ exists, and 

( 10 . 1 ) \{aeaf)^-A\<Ca^^^, 

where C > t), a E (0,1) are constants independent of m, n. If m = 0, we 
have 

(10.2) \{a^afY - A\<Ca^. 

(b) If, in addition, hm|e-/|^oo(o'eO'/)^ = 0, then (7 > 0, a G (0,1) may he chosen 
such that 

(10.3) Ko-eO-/)! < 

10.2. 1-2 and Ising models with period fc x 1. We discuss a special case of pe¬ 
riodic 1-2 models, namely, models with period A; x 1. By the forthcoming Theorem 

10.2, the spectral curves of such models can intersect the unit torus only at real 
points. The conclusions of Theorem 10.1 follow whenever the corresponding char¬ 
acteristic polynomial P{z,w) satishes P(±l, ±1) ^ 0. Similar arguments are valid 
for periodic Ising models, as exemplihed in the forthcoming Example 10.4, which 
illuminates the differences between the assumptions of the current paper and those 
of [29]. 

Let k > 2 and consider a A: x 1 periodic 1-2 model. Let Ha be as in Section 5, 
and let Ha,^,! be the quotient graph of Ha under the weight preserving action of Z^. 
Note that Ha,^,! is a hnite graph that can be embedded in a torus. We can divide 
HA,fc,i into k parts, each of which is bounded by a quadrilateral region enclosing a 
NE/SW edge of the original lattice H, see Figure 10.1. 

There are two vertices of H lying in the Ath such quadrilateral region (see Figure 
5.1), and to these we assign local weights (oji, Cn) and ( 0 * 2 , &i 2 , ^ 2 ), respectively. 
As in Section 5.2, we derive the modihed weighted adjacency matrix K{z,w) of 
Ha,?;,!, with characteristic polynomial P{z,w). 
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Figure 10.1. The decorated graph Ha, 2 , 1 ; the hexagonal lattice H is 
represented by green lines. 


Theorem 10.2. The only possible intersection of the spectral curve {{z, w) : P{z, w) = 
0} with the unit torus = {{z,w) : \z\ = 1, |tc| = 1} is a single real point. 


Proof. Note that a double dimer conhguration is a union of cycles and doubled edges 
with the property that each vertex is incident to exactly two present edges. An 
(unoriented) z-edge is an edge which has weight z or z~^ when oriented. Each 
double-dimer conhguration of the determinant P{z,w) falls into exactly one of the 
following two cases: 

1. it occupies each z-edge exactly once, i.e., each z-edge is in a cycle of the 
double dimer conhguration. 

2. each z-edge is either unoccupied or occupied exactly twice, i.e., each z-edge 
is either absent or a doubled edge in the double dimer conhguration. 

Each k X 1 fundamental domain is comprised of fc 1 x 1 blocks. For 1 < i < u, let 
be the two vertices of the hexagonal lattice lying in the Ah block, see Figure 

10 . 1 . 


Each conhguration in Case 1 consists of a single essential cycle of even length, 
together with some doubled edges. The partition function of conhgurations in Case 
1 is 


1 


Pi — - [AjW -|- Bi] -|- z 


2=1 


2 = 1 


- + B, 

w 
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where 

Ai = — a'f^ai2hi2 — ajia^2^jl — 0*16416^2 + ~ + fli2^i2Cii, 

Bi = 0‘ii0‘i2Ci2 0‘il0‘i2('il “1“ ®ilQlCj2 “1“ ®j2^ilCi2 ®i2Cj]^Cj2- 

Configurations in Case 2 depend on configurations of each 1x1 block, which 
are determined by configurations of boundary edges. Let (respectively, Qp) 
denote the partition function at the ith block when both its x-edges are unoccupied 
(respectively, occupied). Let (respectively, Q^"^) denote the partition function at 
the ith block when its left (respectively, right) x-edge is occupied twice, while the 
right (respectively, left) edge is unoccupied. Then we have 


= 0*10*2 


bilCi2 


+ biibi2 + C*iC*2 — bi2CiiW 


+ biibi2 + C*iC*2 - biiCi2W , 


X I (lil(li2 


Qi — W a^2^iiCii + aiiCi2ai2bii + aiibi2ai2Cii + 6^2^22 [^n — buCn + —j + c^ij , 
= W CL^ibi2Ci2 + 0*2C*iO*i6*2 + 0*26*l0*iC*2 + (^^2 ~ ^i2Ci2 ^ , 

QP = —biibi2CiiCi2 ^ + (0*16*2 + 0*26*1)^ + {OjiiCi2 + 0*2C*i)^ 

+ (o^2^*lC*i + afibi2Ci2 + 0 *i 0 * 26 *lCi 2 + flil(Ii 26 * 2 C*l) ^ , 

where hh = (to — to“^). 

Let X* be the x-edge connecting the tth 1x1 block and the [(i +1) mod /cjth 1x1 
block, and let f* G {0,1, 2} denote the occupation time of x* by a given double dimer 
configuration, which is to say that 

{ 0 if Zi is absent, 

2 if Zi is a doubled edge, 

1 if Zi is in a cycle. 

The partition function in Case 2 is 


E EqI 

tl,...,46(0,2} 1=1 


Let w = e*'^, so that Q‘f‘ > 0. By periodicity, in each configuration, the number 
of (2°^ blocks equals the number of blocks. Therefore in all terms in P 2 where 
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sin 0 appears, which are exactly the terms where and appear, sin 0 has even 
degree. Moreover, given that all the local weights are strictly positive, each term in 
the expansion of P 2 with sin 0 is non-negative. Therefore, 

k k 

p{z, i«)=F,+A=F,+ n + n <2?+ 

i=l i=l 

where F{w) > 0 is the sum of all terms in P 2 in which and appear at some 
1x1 blocks. Let G{z,w) = z[P{z,w) — F{w)]. For given w, G{z,w) is a quadratic 
polynomial in z. Let 

c„=n +F.], c.=n Qf +n of. 

2=1 2=1 2=1 


The roots of G{-,w) are 

(10.4) 

Let tc = and note that 


= 


Ci±y/C'?-4|C'o|2 


2 G. 


Cl - 4|CoP > 4 HiQl^QT] - n + 2AB, 


COS 


v 2 =l 


2 = 1 


Moreover, 


- (A + A + ‘2ABi COS0) 



X {alA2Ci2 + ai2Ciiaiibi2 + ai2&iiaiiQ2 + fenQi c| - 26 ^ 2^2 cos 0 )) 

X {al2^iiCii -|- CLiiCi2ai2bii + (iiibi2(ii2Cii + b^cAbh — 2 biiCii cos 0 — A)) . 


Let ai,bi,Ci > 0. Then Cf ~ 4|CoP > 0 with equality only if w is real. If Gf — 
4|C'oP > 0, then |x^| 7 ^ 1 and G{-,w) has no zeros on the unit circle. Hence G{z,w) 
has no zeros on T x (T \ {±1}). Since there exists at least one pair 9,t E {0,1} with 

P{i-lY, (-ir) = det (- 1 )^) > 0 , 

we have that G{z,w) > 0 on T^, with equality only if w is real. By (10.4), when 
w is real and G^ — dAoP = 0, then is real. Therefore, the only possible 

intersection of P{z,w) with is a single real point. □ 
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The above technique applies also to the spectral curve of the Ising model (not 
necessarily ferromagnetic) with fc x 1 periodicity on the triangular lattice. The details 
are omitted. 

Proposition 10.3. Consider a periodic Ising model on the triangular lattice. To 
each edge e, associate a coupling constant J^, G M. Assume the coupling constants are 
translation-invariant with period k x 1 where k >1. The only possible intersections 
of the spectral curve with the unit torus are real points. 

Example 10.4. Here is an example which explores the generality of the arguments of 
the current paper. We start with an Ising model on the triangular lattice T, with edge 
interactions Je G M that are periodic with period 2x1. Let El he the dual hexagonal 
lattice ofT. The corresponding Fisher graph F is obtained from El by replacing each 
vertex by a triangle. Each triangle edge is assigned weight 1, and a non-triangle edge 
crossing an edge e of T has weight . 

The dimer model on F with the above edge-weights corresponds to the Ising model 
on the triangular lattice. Note that the spins of the Ising model are placed at centres of 
the dodecagons o/F. Two adjacent spins have the same state (respectively, opposite 
states) if and only if the corresponding non-triangle edge of F separating the two 
dodecagons are present (respectively, absent). 

See Figure 10.2 for an illustration, where Oi, 6 i, Ci, 02 , C 2 the edge-weights 



Figure 10.2. The Fisher graph with 2x1 periodic edge weights 

Given the clockwise-odd orientation of Figure 10.2, we can compute the charac¬ 
teristic polynomial P{z,w) = det K{z,w), where K{z,w) is the modified weighted 
adjacency matrix whose rows and columns are indexed by vertices in the 2 x 1 fun¬ 
damental domain. By Corollary 10.3, when Oi, 61 , Ci, 02 , C 2 > 0, the only possible 
intersection of P{z,w) with the unit torus are real. 

The function P{z,w) may be calculated as in Section 5.2, and it may be checked 
that, when (oi, 61 , Ci, 02 , 62 , C 2 ) = (1.1, 0.9, 0, 0, 0.5,0.5), 

(10.5) P(l, 1 )P( 1 , -1)P(-1,1)P(-1, -1) ^ 0 . 
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Let (ai, 6i, 62 , C 2 ) = (1.1,0.9, 0.5, 0.5) and assume that ci, 02 are positive and suffi¬ 
ciently small that (10.5) continues to hold. By Proposition 10.3, the spectral curve 
does not intersect the unit torus. ^45 in [3, 9, 23], the Ising free energy may he ex¬ 
pressed in the form (6.10). However, when (oi, 61 , 62 , C 2 ) = (1.1, 0.9, 0.5,0.5) and 
Cl, 02 are positive and sufficiently small, then neither the high-temperature nor the 
low-temperature condition of [29] is satisfied. 

Moreover, using the technigue of [24] , the sguare of the spin-spin correlation may 
be expressed as the determinant of a block Toeplitz matrix. By applying similar 
technigues as in Sections 8-9, we can obtain the convergence rate of the spin-spin 
correlation of the Ising model of (10.1) whenever the spectral curve does not intersect 
the unit torus. 

10.3. Harnack curve. We present next another sufficient condition for the spectral 
curve P{z,w) = 0 to intersect the unit torus at only real points. The exponential 
convergence rate for two-edge correlation functions follows by Theorem 10.1. 

Harnack curves were studied in [30, 31]. Simply speaking, a Harnack curve is the 
real part of a real algebraic curve A (real zeros of a Laurent polynomial with real 
coefficients) such that the map (10.6) from A to is at most two-to-one. It was 
proved in [22, 23] that the spectral curve of any positive-weight, bi-periodic, planar, 
bipartite dimer model is a Harnack curve. Using the combinatorial results of [10], 
we infer that the spectral curve of any ferromagnetic, bi-periodic, planar Ising model 
is also a Harnack curve. If the Ising model is not ferromagnetic, there may exist a 
concrete counterexample in which the spectral curve is not Harnack, [19]. We give 
here a simple proof that, under certain conditions, the assumption that the spectral 
curve is Harnack implies that its intersections with the unit torus are necessarily 
real. 

Proposition 10.5. Let P{z,w) be a Laurent polynomial taking real values on the 
unit torus. If A := {{z,w) G : P{z,w) = t)} is a Harnack curve, then A can only 
intersect the unit torus at real points. 

Proof. Dehne the logarithmic Gaussian map 'yp \ A ^ CP^ by 



and also Log : H —)■ by 

( 10 . 6 ) 


Log : [z, w) H->• (log 1^1, log |tc|). 


By [30, Lemma 5], H is Harnack, whence the real zeros satisfy H fl = 7 p^(MP^). 
By [30, Lemma 3], ■jf^iM.P^) consists of the singular points of the map Log. 
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Let (z,w) 
that 


and hence 


(e*®,e*'^) E An T^. Since P{z,w) is a Lament polynomial, we have 


dP 

dz 


= 0 , 


dP 

dw 


= 0 , 


'ypiz,w) 


f.dP .dP\ 


Given that P takes real valnes on T^, 7 p(^, w) G Hence H fl C 7 p^(MP^) = 

A n Therefore, any zero of P{z, w) on is real. □ 
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